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Preface 


With  this  thesis,  I  hope  to  help  open  the  book  on 
understanding  the  behavior  and  the  means  of  controlling  per¬ 
iodic  systems.  Many  problems  in  the  field  of  astronautics 
deal  with  periodic  systems  and  as  we  stand  on  t la e  brink  of 
the  exploitation  of  outer  space,  I  feel  we  need  more  under¬ 
standing  of  these  periodic  systems. 

I  would  like  to  extend  my  appreciation  to  my  advisor, 
Dr.  Robert  Calico.  With  his  insight  and  knowledge  I  was  able 
to  overcome  many  problems  encountered  during  the  course  of 
this  study. 

Also,  I  would  like  to  thank  my  wife  Donna.  No  study 
of  this  magnitude  is  without  its  sleepless  nights,  periods  of 
enlightenment  and  times  of  frustration.  She  experienced  these 
times  with  me  and  never  once  displayed  dissatisfaction  over 
the  fact  that  I  was  not  with  her  in  thought  while  at  home. 

Finally,  I  would  like  to  dedicate  this  report  to  my 
father,  who  passed  away  during  its  final  preparation.  It  was 
the  values  he  instilled  within  me  that  made  me  able  to  finish 
this  task. 

G.  Scott  Yeakcl 
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Abs t r  a  c t 


'The  linearized  equations  describing  the  attitude 
motion  oi  two  generic  satellite  cases  are  developed.  in 
both  cases,  the  symmetric  satellite  in  an  elliptical  orbit 
and  an  unsvmmetric  satellite  in  a  circular  orbit,  the  lin¬ 
earized  equations  are  periodic.  Using  Floquct  theory,  the 
stability  oi  the  at  tude  motion  for  several  satellite  design 
is  checked.  For  satellites  with  unacceptable  attitude  stabi¬ 
lity,  an  active  control  scheme  utilizing  modal  control  design 
techniques  is  developed.  With  this  control  scheme,  various 
satellite  test  cases  are  stabilized  and  the  results  arc 
verified  via  a  digital  simulation  of  the  satellite  attitude 
motion. 
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I  MRODUCT  IO\ 
l’ac  ki;  round 

With  the  advent  of  artil'ical  satellites  in  the  late 
1950's,  satellite  designers  have  sought  various  methods  to 
control  the  orientation  of  satellite  bodies.  The  need  for 
proper  attitude,  or  line  of  sight,  control  is  evident  when 
one  looks  at  the  various  missions  that  satellites  perform. 
Whether  a  satellite  is  performing  a  space  surveillance 
mission  for  one  of  the  'super-powers'  or  providing  a  com¬ 
munications  link  for  a  third  world  country,  the  proper 
attitude  must  be  maintaineu  at  all  times  to  assure  the 
satellite  is  accomplishing  its  intended  functions.  A  simple 
example  is  a  reconnaissance  satellite  in  a  low  earth  orbit 
at  an  altitude  of  two  hundred  miles.  If  t lie  line  of  sight 
of  the  cameras  on  board  this  satellite  is  perturbed  by  as 
little  as  five  degrees,  the  camera  will  miss  its  desired 
aimpoint  by  approximately  seventeen  miles.  This  is  totally 
inadequate  performance  for  the  servcillance  of  a  battlefield 
with  high  resolution  cameras. 

To  control  the  attitude  of  satellites,  satellite 
designers  have  used  several  devices  ranging  from  passive 
techniques  using  gravity  gradient  stabilization  to  active 
devices  using  control  moment  gyros,  dual  spin  configurations, 
magnetic  torquers,  or  gyrocompasses.  A  desirable  feature 
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of  gravity  gradient  stabilization  is  that  gravity  torque 
stabilizes  the  satellite  in  lieu  of  a  torque  producing  de¬ 
vice.  Since  no  torque  device,  and  its  necessary  control 
equipment,  is  required,  the  satellite  is  less  costly  to 
build  and  to  launch.  However,  there  are  some  serious  draw¬ 
backs  to  gravity  gradient  stabilization,  which  include 
restricted  satellite  designs  and  low  eccentricity  orbits 
(Kef.  S:  Id'.').  These  two  restrictions  cause  problems  with 
the  physical  design  of  satellites  and  their  intended  mis¬ 
sions.  In  addition  to  these  restrictions,  the  natural  fre¬ 
quency  of  oscillation  of  a  stable  gravity  gradient  satellite 
is  very  low,  and  precise  attitude  control  is  not  attainable 
by  this  means  alone. 

In  the  19b0's,  as  t he  use  of  satellites  increased 
and  the  exploitation  of  space  with  unmanned  vehicles  began, 
this  general  problem  of  attitude  stability  of  satellites 
under  the  influence  of  a  gravity  gradient  was  recognized. 

Two  papers  published  in  the  1900's  dealt  with  the  topic  of 
attitude  stability  of  spinning  satellites  under  the  in¬ 
fluence  of  a  gravity  gradient  torque.  The  first  paper 
examined  the  attitude  stability  of  unsymmctric  satellites 
in  circular  orbits  (Ref.  b:  114).  In  their  study,  Kane 
and  Shippy  examined  satellite  attitude  motion  and  derived 
linearized  equations  with  periodic  coefficients.  Since 
the  linearized  equations  had  periodic  coefficients,  their 
asymptotic  stability  was  determined  by  using  Tloquet  theory. 
In  fact,  the  stability  of  this  class  of  satellites  was 


determined  to  be  a  function  of  moment  of  inertia  ratios 
and  spin  rates  (Ref.  0:  1  1 4  -  1  1  ‘J  J  .  In  the  second  study, 

Kane  and  Barba  examined  the  stability  of  symmetric  satellites 
in  elliptical  orbits.  Like  the  first  study,  this  study 
provided  a  periodic  linear  system  model  of  the  satellite 
attitude  motion.  Once  again,  i-ioquet  theory  was  used 
to  determine  the  stability  of  this  class  of  satellites. 

During  this  study,  stability  was  found  to  be  dependent  on 
the  spin  rate  and  moment  of  inertia  ratios  like  the  first 
study,  and  on  the  orbital  eccentricity  (Ref.  7:  402-405). 

These  two  studies  were  concerned  exclusively  with  the  passive 
stability  of  the  satellite  attitude;  there  was  no  attempt 
to  control  those  satellites  because  of  the  non-constant 
coefficients  of  the  linearized  systems. 

The  results  of  these  two  studies  provided  valuable 
information  for  satellite  engineers  in  the  llKul's  and  1070's, 
liuwever,  as  we  become  more  and  more  dependent  on  satellites 
and  find  new  uses  for  satellites,  we  may  no  longer  be  able 
to  satisfy  the  requirements  for  passive  stabilization  or  we 
may  simply  want  to  increase  the  low  damping  of  tire  natural 
low  frequency  oscillations.  In  fact,  with  the  advent  of  the 
Space  Transportation  System,  many  space  planners  arc  devclopin 
concepts  of  large  space  structures  for  use  in  the  future.  One 
of  the  general  mission  requirements  of  these  large  space  struc 
turcs  is  precise  vehicle  line  of  sight  stability  (Ref.  4:  52). 
With  this  requirement,  considering  that  low  orbital  ccccntrici 
tics,  spin  rates,  and  vehicle  design  requirements  need  to 


bo  overcome,  an  active  control  scheme  for  these  satellites 
must  be  developed.  However,  a  problem  arises  in  the  devel¬ 
opment  of  an  active  controller  for  these  satellites.  The 
linear  systems  used  to  model  the  attitude  dynamics  of  both 
the  unsymmetric  satellite  in  a  circular  orbit  and  the  symme¬ 
tric  satellite  in  an  elliptical  orbit  are  nonautonomous ,  or 
specifically  periodic,  systems.  As  such,  these  systems  do  not 
allow  the  use  of  classical  linear  controller  design  techniques. 

The  control  of  periodic  systems  is  a  topic  which  is 
in  need  of  further  research.  Many  references  can  be  found 
on  the  control  of  constant  coefficient  linear  systems.  How¬ 
ever,  ver;  little  material  can  be  found  on  periodic  systems 
as  well  as  the  control  of  such  systems.  In  his  study  in 
1980,  Shelton  (Ref.  141  found  modal  control  techniques 
feasible  to  control  the  orbit  of  a  satellite  in  orbit  about 
the  third  libration  point,  where  the  orbit  is  a  time- 
periodic  system.  Therefore,  since  this  type  of  control 
methodology  was  used  to  control  a  time-periodic  system,  it 
is  used  in  this  study  for  the  control  of  satellite  attitude 
motion . 

Modal  control  is  a  modern  control  technique  which  at¬ 
tempts  to  independently  control  systems  modes  by  means  of 
providing  state  feedback  which  is  proportional  to  the  eigen¬ 
vectors  of  the  dual  basis  for  the  system.  Using  direct 
feedback  of  the  modal  variables,  tire  closed-loop  eigenvalues 
of  a  constant  coefficient  system  can  be  placed  at  a  desired 
location  using  pole  placement  techniques.  Also,  when  a  con- 


stant  coefficient  system  is  controlled,  the  time  response  of 
the  system  can  be  controlled  directly  through  the  gain  selection. 
In  this  study,  modal  control  techniques  are  used  to  control 
the  location,  in  the  complex  plane,  of  the  characteristic  ex¬ 
ponents.  These  characteristic  exponents,  of  both  the  open 
and  closed- loop  systems,  arc  found  using  Floquet  theory. 

Problem  and  Scope 

The  problem  addressed  by  this  study  is  the  control 
of  linear  systems  with  periodic  coefficients.  Shelton 
developed  a  controller  for  a  particular  periodic  system  in 
his  research.  However,  the  system  lie  addressed  had  one 
unstable  characteristic  exponent.  To  handle  this  system, 
lie  developed  a  simplified  controller  to  control  only  one 
mode.  In  this  study,  the  general  form  of  a  controller  cap¬ 
able  of  controlling  any  number  of  modes  is  developed.  This 
general  controller  is  implemented  in  several  example  cases 
and  the  results  are  verified  by  the  use  of  a  digital  simul¬ 
ation  of  the  linearized  closed-loop  system. 
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CHAPTER  II 


PROBLEM  ANALYSIS 

The  major  thrust  of  this  study  is  to  develop  a 
control  scheme  for  the  use  with  linear  periodic  systems.  As 
mentioned  previously,  the  two  periodic  systems  of  interest 
are  the  symmetric  satellite  in  an  elliptical  orbit  and  the 
unsymmetric  satellite  in  a  circular  orbit.  In  this  chapter, 
the  attitude  dynamics  of  these  two  satellite  cases  arc 
examined.  With  these  results,  the  linearised  equations  for 
the  satellite  attitude  motion  arc  developed. 

Satellite  Attitude  Kinematics 

The  first  step  in  developing  the  equations  of  motion 
for  the  attitude  dynamics  of  a  satellite  is  the  derivation 
of  the  kinematics  of  the  problem.  This  step  involves  the 
establishment  of  coordinate  frames ,  angular  displacements , 
and  angular  rates.  As  depicted  in  figure  1,  the  basic 
orbital  layout  consists  of  a  satellite  S  orbiting  body  F, 
with  an  assumed  inertial  frame  at  F.  Additionally,  n  and  v, 
which  represent  the  range  and  true  anomaly  of  the  satellite, 
are  the  polar  coordinates  of  5  with  respect  to  F. 

The  angular  orientation  of  satellite  S  in  inertial 
space  is  found  through  a  series  of  rotations.  The  first 
rotation  to  Frame  1  is  a  rotation  equal  to  the  true 
anomaly  about  the  unit  vector.  Subsequent  rotations  to 
a  body-fixed  frame  arc:  a  rotation  of  0  about  the  i  axis, 

i  ) 


b 


a  rotation  of  d 5 about  the  j  axis,  and  a  rotation  of  ^  about 
the  k  axis.  Those  angles,  0,0,  and  9  ,  represent  the 

3  12  3 

pitch,  roil,  and  yaw  of  the  satellite  body. 

Kith  these  angular  rotations,  the  total  angular  ve¬ 
locity  and  acceleration  of  S  with  respect  to  inertial  space 
can  be  determined.  by  using  the  angular  velocity  chain 


rule , 


(  I 

-43  -32 

03  '  +  u) 


(1) 
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•  •  ■  ..  ’  • 

0  cosO-O  0  sinO  -vcos  0  sin  0  +  0  0  +  v 
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With  those  two  vectors,  the  angular  velocity  atul  angular 
acceleration  of  S  with  respect  to  inertial  space,  the  equa¬ 
tions  tor  the  satellite  attitude  dynamics  can  lie  derived. 


Symmetric  Satellite  in  an  l.lliptical  Orbit 
(Ref.  7:  402-404) 

The  development  of  the  equations  of  motion  for  tiic 
attitude  of  a  symmetric  satellite  in  an  elliptical  orbit 
perturbed  by  a  gravity  gradient  torque  is  given  by  Kane 
and  Barba  (1900).  The  following  derivation  is  adopted  from 
this  reference. 

The  first  step  in  this  development  is  the  investi¬ 
gation  of  the  orbital  motion  of  a  satellite  in  the  common 
two-body  orbit,  figure  1.  In  additon  to  the  parameters 
specified  in  Figure  1,  the  orbital  eccentricity  is  noted 
as  e.  and  a  is  the  semi-major  axis  of  the  orbit.  In  the 
polar  coordinates  A  and  v,  the  differential  equations 
describing  the  two-body  orbit  are, 

A  -  Av 2  +  11  =  0  ( 6a ) 

A7 

and 


A2v  =  sja.  y(l-c2)‘! 


(bb) 
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where  p  is  the  Gravitational  Parameter.  The  two-body  satellite 
motion  is  periodic,  and  the  period  may  be  found  from  Kepler's 
Second  Law  (Ref.  1:  55): 


*  tv  / 

With  this  result,  liquation  (0)  can  be  rewritten 


and 


/i2v 


(7) 


(8a) 


(Sb) 


Substituting  (Sb)  into  (Sa) ,  one  obtains  a  differential 
equation  for  the  range  of  a  satellite,  in  an  elliptical  orbit, 
from  the  central  body: 


K  -  H 


—  (W2)  V 

1  A2 


(9) 


This  differential  equation  can  be  nond imensional i zed  by 
making  the  following  variable  transformations: 


where  t  and  g  represents  nond imens iona 1  time  and  distance, 
respectively.  With  these  nondimensional  variables,  liquation 
(9)  becomes 


I  T 


(c2-D 

b3 


0 


(10) 


where  the  prime  indicates  derivatives  with  respect  to  x.  The 
resulting  liquation  (10)  is  a  second  order  differential  equation. 
As  such,  it  requires  two  initial  conditions:  Assuming  that 
the  satellite  begins  at  perigee. 


^perigee 


a(l- c) 


the  initial  conditions  tor  equation  (10)  arc 


' ( u )  =  1-c 
' t  o )  =  0 


At  this  point,  the  equations  of  motion  tor  the  satel¬ 
lite  attitude  can  be  derived  using  liuler's  Moment  liquation: 

{Ml  =  L  1  J  (q  c.1111}  +  {5s,i)ri]  fw"*1}  (11) 

3  3  OT  3  3  3  3 

where  [1J  is  the  mass  moment  of  inertia  tensor 

[uij  is  the  shew  matrix  associated  with  o  I  *  ,  and 
[M]  is  the  external  moment  vector. 

Assuming  that  the  axes  are  defined  in  such  a  way  that  frame 
3  lies  along  t lie  principal  axes  of  the  satellite  body,  the 
inertia  tensor  takes  the  form 


1 


i  i 


J 


[0] 


0] 


1 


(12) 
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Additionally,  the  skew  symmetric  matrix  associated  with 
angular  velocity  has  the  form 


3  i 


0 


-  id 


(13) 


,  and  are  the  individual  components  of  m 

12  3 


where 

Substituting  liquations  (3),  (3),  (12),  and  (13)  into  1 1  q u a 
tioti  (11)  produces  moment  equations  in  the  form: 
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{M}  = 

liquating  liquations  (14)  and  (15),  recalling  that  for  a 

symmetric  satellite,  I  =1  ,  one  finds 

11  2  2 


I  oj  '  +  (I  -I  )u  a  =  0 

111  3  3  2  2  2  3 

I  CO  '  +  (  I  -I  )  CO  CJ  =  2) 

112  11  3  3  13 


V* 


cos©  sin6 
2  2 


03  '  "0 


Finally,  by  defining  an  inertia  parameter  k  such  that 

i 


I  -  I 

3  3  ; 


1  l 


liquation  (lb)  becomes 


(15) 


(16) 


o) '  +  k  co  co  =0 

1  12  3 


-  k 


■)  + 


1  1 


k  cosd  sine 


=  0 


co’  =  0 

3 


(1") 


The  three  nonlinear  differential  equations  given  by 
liquation  (17)  can  be  combined  with  liquation  (5)  to  obtain 
three  second  order  equations  in  the  variables  0,0,  and  0  . 

1  2  3 


By  do  fining  the  state  vector  to  he 


xU) 


VTl 

0  i  T  ) 


o;u-) 


o;  (x) 

and  noting  the  reference  state. 


x  =  0 

c. 

a  linear  system  of  the  form 
x'  (.t)  =  A  (. T  )  X  ( T  ) 


can  be  derived  in  the  neighborhood  of  tiie  reference  state. 
To  obtain  this  linear  system,  the  original  system 


x  '  (.T  )  =  l'[x  U  )  ] 

is  linearized  about  x  such,  tiiat 

o 

Aij  =  X  =  X  1,J  =  1,  ....  4 

J  o 

The  resulting  coefficient  matrix,  A,  is  a  function  of  the 
inertia  parameter,  k  ,  and  the  variables,  \> '  ,  v",  and  O'. 

1  3 

To  produce  a  usable  linear  system,  the  dependence  on  these 

variables  must  be  eliminated.  The  dependence  on  v'  and  v" 

can  be  eliminated  through  the  use  of  equation  (.Sb)  .  From 

Hquat ion  (.Sh)  , 


v  ' 


C 1  *  c 2 ) '" 


(.IS) 


and  by  nond  imens  ional  i  a  ing  with  and  x,  one  obtains 


which  by  d i  fferent iat ion  with  respect  to  i  yields, 

,  I 

V  =  :  ;•  (20) 

,  3 

Replacing  v'  and  v”  with  liquations  (ID)  and  (20)  eliminates 
the  dependence  on  these  variables.  To  eliminate  the  de¬ 
pendence  on  one  recognizes  from  liquation  (10) 

3 

U)  '  =  o 

3 

implying  that 

m  =  v’  +  O’  =  constant 

3  3 

The  dependence  on  9'  is  eliminated  by  specifying  this  constant 
However,  the  constant  represents  the  spin  rate  of  the  satel 
li.tc  and  the  true  anomally  rate.  Since  one  is  concerned  only 
with  the  spin  rate  of  the  satellite,  the  constant  is  defined 
as  (a  +  i)  and 

v  '  +  u  '  =  (u  +  1 )  (21) 

Substituting  liquations  (19)  through  (21)  into  the  linearized 
coefficient  matrix,  a  linear  equation  is  found  in  the  form 
o  f 

X  ’  (x  )  =  A  (T  j  X  (T  ) 

where  A  ( t )  is  explicitly  a  function  of  and  implicitly  a 
function  of  t.  The  individual  elements  of  this  coefficient 
matrix  for  the  linearized  system  are  presented  in  Table  1. 

Upon  examining  Table  I,  one  notes  the  coefficient 
matrix  is  a  function  of  the  orbital  eccentricity,  the  spin 
parameter,  the  inertia  parameter,  2  and  t.  '  .  Since  g  repre- 


seats  the  nond imcns ional  distance'  from  the  satellite  to 
body  F,  it  is  periodic  with  the  same  period  as  the  orbital 
period.  Thus, 

'  lit  =  '.l  T  +  ) 

where  T  =  1 ::  and  since  A  is  a  function  of  4 

A(t)  =  A  ( •[  +  7)  (22) 

Thus,  the  equations  for  the  attitude  motion  of  the  symmetric 
satellite  are  periodic  and  this  periodicity  is  due  to  the 
orbital  motion  of  the  satellite. 

Unsymme t  r  i  c a  1  Sa t  e 1 1  i t  e  in  a  Circular  Orbit 
(Ref.  6:  114-117) 

The  development  of  the  equations  describing  the 
attitude  motion  for  the  unsvmmetric  satellite  in  a  circu¬ 
lar  orbit  perturbed  by  gravity  gradient  torque  arc  given 
by  Kane  and  Shippy  (1905),  from  which  the  following  de¬ 
rivation  is  adopted. 

Unlike  the  symmetric  satellite,  this  case  does  not 
require  analyzing  the  orbital  motion.  Therefore,  the  first 
step  in  the  derivation  is  to  employ  Ruler's  Moment  liquation 
(11).  In  this  case,  since  the  satellite  is  not  symmetric 
the  body  fixed  axes  are  used  as  the  reference  frame.  Using 
coordinate  t  rails  forma  t  ions ,  liquation  (5)  can  be  written  as 


10 


10  coso  -vcoso  sino  )cose  + 

12  12  3 

(0  +vsinO  )sino 

2  1  3 

{</'  i)  =  (.0  +vsinO  )cosO  -  (6  coso  -  (23 

2  13  12 

osimj  COSO  )sinO 

2  1  3 

Q  +0  sine  + vcosO  cosO 

3  12  2  1  J  4 

Since  the  attitude  ol'  the  satellite  in  this  case,  like  the 

symmetric  satellite,  is  a  function  of  0  and  0  ,  only  small 

1  2 

angles  are  considered.  Considering  only  small  angles,  liqua¬ 
tion  (25]  can  be  linearized  about  these  small  angles  by 


using  the  following  approximations: 


co so  =  1 

1 

eosO  -  1 

2 


s i no  -  0 


Linearizing  liquation  (25),  one  obtains 


(G  -vO  )coso  +  (6  +v0  )sino 

1  2  3  2  1  3 

{ a>  'Z1]  =  (0  +v0  )cos0  -  (6  “v G  )sinO 

2  3  3  1  2  3 

0  +  V 

3 

Since  the  orliit  is  circular,  v  is  constant  and  v  is  zero, 
for  all  time.  Denoting  v  as  i,  the  angular  acceleration 


is  found  from  liquation  (25): 


(0  -AO  H'Osq  -(0  -AO  )sino  0  + 

2  3  3 


l 

2  3  1 

(0 

+ A  ■*  isino  +(6 

2 

1  3 

,  »  k  I  i . 

•  • 

to  ' 
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moments  acting 

on 

5 

( 1  inear 

ized  by  using 

approximat ions 

o  f 

liqua  t  ion 

(24))  are: 
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3  J 


(2: 


With  tlic  Huler's  Moment  liquation  (11)  along  with  liquations 
(25)  through  (27)  and  the  following  inertia  parameters: 


J 


33 _ l  2 
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1  -  I 

k,  =  --1 - -- 
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I 


22 


1  1 


122 


1  1 
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3  3 


The  equations  describing  the  attitude  motion  are  given  by: 


sinG  {4k  A20  +0  -0  0  (1+k  )+A0  (1-k  )+A0  0 

3  1  2  2  13  1  1  1  3  2 


(1+k  ) } +cos0  (0  -A0  (1+k  ) +0  0  (1+k  )+A0  0 

1  3  11  12  3  1  13 


(2  3; 


(1+k  }+k  A  ~  0  } 
i  l  i 


si nu  (-0  +A0  (1  +  k  ) - 0  0  (1-k  )-A0  0  (1-k  ) 

3  1  2  2  2  3  2  1  3  2 


+  k  A 2  0  }  +cos0  (-4  A^k  0  -0  0  (1-k  )+0  +A0 

2  1  3  2  2  1  3  2  2  1 


(231 


(1+k  ) +A0  0  (1-k  )}  = 

2  2  3  2 


0 


U  -3A:'k  si  no  cosO  =  0  ( 

3  3  3  3 


Defining  a  parameter,  p;  such  that 

(k  +  k  ) 

->  _  -  A  2  1  ? 

-i+TTk— 

1 


liquation  (28c)  becomes 


0  +  p 2  s  i no  c o s 0  =  0  ( 

3  3  3 


This  second-order  differential  equation  for  0  (t)  contains 

3 

a  first  integral 


•  n -  ,  i, 

0  -  r  L  (cos20  - 1 )  +  0  - ]  '  ( 

3  ~  -  3 

where 


0  =  0  3  ( t  =  0 ) 


Defining  a  new  independent  variable,  x,  such  that 


t  =  20  ( 

3 

liquation  (50)  can  be  rewritten 

0=  +  [^-(eosT-l)+02]"  ( 

3 

A t  this  point,  the  linearized  equations  describing  the 
attitude  motion  can  be  obtained.  By  defining  the  state 


vector : 


and  a  function  ft,  such  that 


3  =  (34) 

a  linear  system  of  the  form  of 
X  ’  (  X  )  =  A  (  X  )  X  1  T  j 

is  found,  flic  differential  equation  for  x  (t)  and  x  (x) 

1  2 

are  found  by  using  liquation  (31),  (33),  and  (34): 

x|  (x)  =  X  s  (x) 
and  likewise, 

X’ (X)  =  4  X  (X) 

’  -  I* 

The  differential  equations  for  x  (x)  and  x  (x)  arc  found 

3  4 

from  linear  combinations  of  liquations  (28a)  and  (28b).  The 
resulting  coefficient  matrix  is  given  in  Table  II. 

The  function  ft  can  be  found  from  liquations  (2Sd), 

(52),  and  (34)  such  that 

k  J  n*1 

-  j)--y  (COST-1)  +  —  (5  5) 

1  2J  A2J 

IV i t h  this  relation,  all  elements  of  Table  II  are  functions 
of  the  independent  variable  x. 


Table  11 


Non -Zero  Coefficient  Matrix  Clements 
for  the  Unsymmetric  Satellite 

A  ( 1 , 5 )  =  | 

A(2,4)  -  V 

A  ( 3 . 1 )  =  -.5  +  [k  (l-cost)-k  ( 1  +  e  o  s  r  )  ] 

‘*2  1 

A ( 3 , 2 )  =  -iiilil  (kj+k2)  sin  t 

A ( 3 , 3 )  =  ^4—  (k  +k  )  sin  t 
4  1  2 

A  (3 , 4 )  =  +  44*  [k  (1-cost)  -  k  ( 1  +  cost  )  ] 

^  L\  Z  1 

A (4 , 1 )  =  - 11411  (k  +k  )  sin  t 
4  1  2 

A  ( 4 , 2 )  =  -.5  +  --4-  [k  (1+cost)  -  k  (.1- cost  4] 

2  1 

A  ( 4 , 3 )  =  -  [k  (1+cosx)  -  k  (1-cost)] 

H  2  1 

((3  +  1) 


A(4 , 4) 


(k  +k  )  sin  t 


By  examining  Table  il,  it  is  apparent  that  all  elements  are 
periodic  with  period  t ,  such  that 

A  (  r  )  =  A  (  t  +  T ) 


w  li  e  r  e  T  =  -  :r . 

hxamining  the  elements  of  the  coefficient  matrix  in 
Table  11,  one  notes  the  dependence  of  A(x)  on  the  three 
dimensionless  parameters  k  ,  k  ,  and  8 .  The  parameters  k 
and  k  arc  specified  by  the  size,  shape,  and  mass  distribu¬ 
tion  of  3 .  The  parameter  B,  however,  is  dependent  on  the 
independent  variable  t  and  on  C2/A%  which  is  a  measure  of 
the  rotational  rate  of  the  satellite.  The  value  of  O'/A2  can 
be  found  for  a  particular  satellite  configuration  and  spin 
rate.  Defining  a  spin  parameter,  a,  and  a  parameter,  k , 
such  that 


a 


(5b) 


where  0  is  t lie  average  value  of  0  during  one  revolution  of 

3  3 

S ,  the  value  of  b2/A“  can  be  found  for  any  value  of  a.  Re¬ 
writing  liquation  (52)  such  that 


i)  =  + 
3  — 


k2  o2 


(COST-1) 


(57) 


and  by  using  the  method  of  separation  of  variables  and  the 
ident i ty , 

■z  (cost-1)  =  -sin2T 


liquation  (57)  can  be  rewritten 


+  0  Jt  = - 3 - j-  (IS) 

3  (1- fc2sin2T)  - 

Defining  T  as  the  time  required  for  0  to  change  by  an  amount 

3 

of  277  radians,  liquation  (58)  can  be  written 


de 


(1- fc2sin?T) ' 


(59) 


'  o  o 

where  the  left  hand  side  of  liquation  (59)  is  in  the  form  of 
complete  elliptic  integral  of  the  first  kind.  Performing 
the  integrations,  one  finds 


4fl(fc)  =  +  TO 


(40) 


Since  0 „  is  an  average  value  during  one  revolution, 

5  27T  Tt  0 


and 


a  = 


jjrrcf:) 


or 


A 


25  (t)^ 


(41) 


Thus,  given  a  desired  ot,  the  parameter  0/A  can  be  found  from 
liquation  (41)  by  finding  k  and  then  5(1).  Substituting 
Equations  (41)  and  ( 2 S d )  into  Equation  (50),  one  obtains 


TT 

2a 


-5(k^k2)  h 

1  +  k  k 
1  2 


5(1) k  = 


+ 


(42) 


which  yields  a  transcendental  equation  in  terms  of  k.  Solv¬ 
ing  this  transcendental  equation  for  liquation  (41)  and 

(55)  can  he  used  to  solve  for  file  function  p  to  complete  the 
solution  of  the  A ( t )  matrix  calculations  in  Table  II. 

This  completes  the  analysis  of  the  unsymmctric  satel 
lite  in  a  circular  orbit  around  body  F.  A  linear  periodic 
system  of  equations  was  found  with  a  period  of  It.  Unlike 
the  first  case,  where  the  periodicity  is  due  to  the  ellip¬ 
tical  orbital  motion,  the  periodicity  for  this  case  is  due 
to  the  rotation  of  the  satellite  about  its  body  fixed  k 


CHAPTER  III 


Fl.OqUliT  THHORY 
Periodic  Systems 

In  Chapter  II  oL'  this  study,  the  linearized  equations 
describing  the  attitude  motion  for  the  two  satellite  cases 
are  derived.  For  both  cases,  the  equations  of  motion  are 
found  to  be  linear  of  the  form, 

x  ’  (T)  =  A  ( t ) x  (x)  (45a) 

wliere  the  coefficient  matrix,  Ait)  is  periodic  such  tliat 

A (t  )  =  ACT  +  T)  (43b) 

with  T  being  defined  as  the  period  of  the  linear  system. 

For  periodic  systems  of  this  type  the  stability  can 
be  determined  through  tiie  use  of  Floquet  theory.  The  essence 
of  Floquet  theory  allows  equations  in  the  form  of  liquation 
(43)  to  be  theoretically  reduced  to  the  case  of  constant 
coefficient  equations  (Ref.  5:  oO) . 

Floquet  Theory 

In  general,  Floquet  theory  states  that  the  fundamental 
matrix  for  ordinary  differential  equations  in  the  form  of 
Lquution  (43)  is  given  by 

4(t)  =  P  (t  )  eFx  (44) 

where  F  is  a  constant  matrix  and  P(t)  is  a  periodic  matrix 
with  the  same  period,  T,  as  the  coefficient  matrix,  A(t) 


(Ref.  12:  2 0 5 )  .  This  fundmental  matrix  is  similar  to  tlie 
fundamental  matrix 

t^(x)  =  (4  5) 

where  F(t)  =  1'(t  +  T)  and  J  is  the  Jordan  form  matrix 
associated  with  the  matrix  F.  Tor  the  case  where  r  has  n 
distinct  eigenvalues  A^,  tlie  Jordan  form  matrix,  J,  is 
diagonal.  The  n  linearly  independent  solution  vectors 
making  up  '!'(t)  in  this  case  have  the  form 

=  fi(T)e>'lT  (46) 

where  the  f.(x)  are  the  column  vectors  of  I-'(t).  The 
eigenvalues,  A ^ ,  of  the  matrix  F ,  which  are  termed  the 
characteristic  exponents,  determine  the  asymptotic  stabil¬ 
ity  of  the  linear  system.  Specifically,  if  all  the  charac¬ 
teristic  exponents  have  negative  real  parts,  the  system  de¬ 
scribed  by  Equation  (43)  is  stable.  If  the  eigenvalues  of  F 
are  not  distinct,  the  Jordan  form  matrix  is  block  diagonal 
and,  like  the  distinct  case,  stability  is  assured  if  all  the 
eigenvalues  have  negative  real  parts.  Therefore,  the  system 
stability  can  be  determined  from  the  A.'s.  It  remains  to  be 
shown  how  these  characteristic  exponents  can  be  computed. 

In  order  to  compute  the  characteristic  exponents  and 
the  periodic  matrix  P(t),  we  consider  the  following.  The 
principal  fundamental  matrix,  or  state  transition  matrix 
over  one  period,  <1>(T,0),  referred  to  as  the  monodromy  matrix, 
is  determined  from  the  solution  to  the  general  matrix 
differential  equation: 


[0(t)  ]  =  A(t)Ht) 


(47) 


U 

Jt 

with  the  initial  conditions 
0(0)  =  [I] 

The  monodromy  matrix  for  a  linear  system  of  the  form  of 
liquation  (45)  can  be  found  by  direct  numerical  integration 
of  liquation  (47).  With  the  monodromy  matrix,  a  general 
solution  to  liquation  (45)  at  c  =  T  is  given  by 

x(T)  =  1 ( T , 0 ) x ( 0 )  (43) 

where  x(0)  is  the  initial  conditions  vector.  From  Equation 
(46),  the  solution  at  t  =  T  corresponding  to  a  specific 
value  of  A^  is: 

x  ■  ^ 

x i ( T )  =  e  1‘fiCT)  (49) 

with  the  initial  condition  vector  given  by: 

xt(0)  =  f^O)  (50) 

Substituting  Equations  (49)  and  (50)  into  Equation  (4S) 
y  i  e  1  d  s  : 

XiT- 

<P( T,0)fi(0)  =  e  fi(T)  (51) 

where  due  to  the  periodic  nature  of  f.(t), 
f i (T)  =  fi(0) 


Rearranging  Equation  (51), 


p  {  B  - 1  p  }  v  =  0 

j  where  p  is  an  eigenvalue  and  v  is  an  associated  eigenvector 

|  X  i'!' 

!  ot  matrix  R.  Thereiore,  bv  analogv,  e  is  an  eigenvalue  and 

j  •  - 

I  f . (0)  is  an  eigenvector  of  the  monodromv  matrix.  Thus,  by 

|  1 

j  evaluating  the  eigenstructurc  ol  the  monodroiny  matrix,  the 

I  characteristic  multipliers,  p.,  defined  by, 

Ai7 

P I  =  e  (in; 

are  found.  With  these  n  characteristic  multipliers,  t lie 
j  characteristic  exponents  are  found  by, 


Using  these  character i st  ic  exponents,  the  Jordan  form  matrix, 
J,  of  liquation  (45)  is  found  such  that  (Ref.  12:  2o7): 


Ia  addition  to  find  ini;  the  characteristic  exponents 
and  the  Jordan  l'orm  matrix,  for  a  complete  solution  in  the 
form  of  liquation  (45),  the  F(x)  matrix  must  be  found  for  all 
t.  Since  F  (r)  is  periodic  with  period  T,  finding  F(x), 
xc(0,T),  is  equivalent  to  finding  F(x)  over  all  t.  Defining 
a  matrix  •;  ( x  )  such  that: 


a  fundamental  matrix  is  formed.  Ana  1 gous  to  this  matrix, 
we  can  form  a  square  n  x  n  matrix  F(t)  sucli  that: 


i;(t)  =  [f  (x)  :  f  ft)  :  •  •  •  :  r  ft)]  (57) 

12 

Substituting  tiie  matrices  F(x)  and  J  into  liquation  (45),  the 
fundamental  matrix,  as  given  by  Equation  (56),  is  found. 

Being  a  fundamental  matrix,  it  follows  the  matrix  differen¬ 
tial  equation  given  by  Equation  (47)  such  that: 

=  F  ( x )  e  J  L 

and 

qVn(x))=  F'(t)eJl  +  F'(T)JeJl  =  A(t)F(t)cJt  (5S) 

Jr 


Multiplying  by  tlic  inverse  of  the  matrix  c 
Equation  (58)  becomes 


,  and  simplifying, 


!■'  Mr)  =  A (t  )  1;  (i  )  -  1;  (t  )  J  (  59) 

where  the  initial  condition  mat  r  ix  ,  f  (  U  )  ,  is  given  by  the 
partitioned  matrix, 

1  C 0)  =  :  ( to  :  f  ( o) :  *  •  *  :rn(in]  (oo) 

Thus,  F  (r )  can  be  found  by  integrating  the  matrix  differen¬ 
tial  liquation  (59),  over  one  period,  with  the  initial  con¬ 
ditions  given  by  liquation  (00)  . 

Kith  these  evaluation  of  the  Jordan  form  matrix,  J, 
and  the  eigenvector  matrix,  F(t),  all  system  information  is 
found  to  complete  a  Fioquet  solvit  ion  or  the  linear  periodic 
system  of  liquation  (43).  In  the  following  chapter,  these 
results  are  used  to  develop  a  control  law  to  provide  active 
control  for  these  periodic  systems. 
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CONTROL  Til  LORY 

In  his  thesis  on  the  control  of  a  satellite  orbit, 
W.L.  Shelton  found  that  modal  control  techniques  provided  a 
feasible  control  scheme  for  time-periodic  systems  (Ref.  14) 
Basically,  modal  control  techniques  allow  the  control  engi¬ 
neer  to  decouple  the  various  system  modes  and  control  them 
individually.  In  essense,  in  a  completely  controllable  con 
stant  coefficient  system,  the  poles  of  the  system  can  be 
placed  where  desired  using  full  state  feedback.  With  these 
pole  placement  techniques,  the  characteristics  of  each  mode 
can  be  specified. 

liven  though  any  number  of  modes  can  be  controlled 
using  modal  control  techniques,  Shelton's  system  had  only 
one  mode  requiring  control.  Thus,  the  controller  developed 
was  simplified  to  handle  one  mode,  which  is  not  the  general 
case.  In  general,  one  would  like  to  control  any  number  of 
modes  simultaneously.  As  shown  in  Chapter  V  of  this  study, 
the  two  example  satellite  cases  at  times  exhibit  more  than 
one  unstable  characteristic  exponent,  or  mode.  Therefore, 
a  generalized  controller  capable  of  controlling  multiple 
modes  is  required.  In  this  chapter,  t  he  design  and  imple¬ 
mentation  of  this  generalized  controller  is  presented. 

Modal  Control  Theory 

The  basic  theory  behind  modal  control  is  that  any 
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system  can  be  decoupled  into  individual  modes.  These  indi¬ 
vidual  modes  are  determined  by  the  characteristics  of  system 
poles.  Therefore,  if  the  system  is  fourth  order,  there 
are  four  modes.  To  decouple  a  system  into  its  individual 
modes  one  needs  to  consider  the  mathematical  properties  of 
similarity  matrices.  Matrices  J  and  A  arc  said  to  be  similar 
i  f 

J  =  f ' 1 Af  (01) 

where  J,  I:,  and  A  are  constant  n-by-n  matrices  with  f  being 
nonsingular.  Since  J  and  A  have  the  same  eigenvalues,  this 
equation  can  be  used  for  pole  placement  (Ref.  9:  61). 
furthermore,  d  can  be  chosen  in  its  simplest  form,  the 
Jordan  canonical  form.  for  n  distinct  real  eigenvalues  of 
A,  the  Jordan  form  matrix  has  the  structure, 

,o 

i 

o 

J  =  •  (62) 


where  the  ct 1 s  represent  the  eigenvalues.  The  f  matrix 
required  for  this  similarity  transformation,  consists  of  a 
partitioned  matrix  of  the  n  eigenvectors  of  A: 


where  IT  is  the  eigenvector  associated  with  a • .  Thus,  with 
the  T  matrix  given  by  liquation  ,  the  matrix  J  is  i'ound 

to  Ik-  similar  to  the  matrix  A.  Applying  this  concept  ,  a 
constant  coefficient  system  can  be  decoupled.  Given  the 
system , 

x  1 1 )  =  Ax  It)  (64) 

the  n  modes  can  be  decoupled  by  defining  a  new  set  of 
coordinates  ,  n  (t)  , 

x  ( t )  =  It,  (t)  (0  5) 

Taking  the  derivative  of  liquation  (.65)  with  respect  to  the 
independent  variable,  t,  one  finds 

x(t)  =  l'n(t)  (06) 

Substituting  liquations  (65)  and  (661  into  the  original 
system  equation, 

i;n(t)  =  ATn(t) 

or 

rut)  =  i;  *  *  Al:r,  1 1 )  (67) 

Thus ,  if  T  is  made  up  of  the  n  eigenvectors  of  A, 

r;  |  t )  =  d  n  ( t )  ( 6  S ) 

This  resulting  system,  expressed  in  modal  variable,  consists 
of  a  set  of  n  uncoupled  differential  equations: 


r'l  (t)  -  o  n  (t) 

1  1  1 

r'l  (t)  =  0  n  (t) 

222  (0‘J) 

Since  the  eigenvalues  of  J  and  A  are  t lie  same,  the  stability 
characteristics  and  responses  of  the  individual  modes  lias  not 
changed.  This  same  technique  of  decoupling  the  system  modes 
can  be  applied  to  linear  periodic  systems  in  the  form  of 
liquation  (45),  using  the  results  of  Chapter  III.  However, 
since  control  of  the  system  is  desired,  a  control  term  must 
be  added  to  the  state  equations. 

Applying  control  to  the  system,  the  state  equations 
in  the  form  of  liquation  (45),  are  augmented  with  a  control 
term, 

X'(t)  =  Aft  )  X  (t )  +  b(T)u(T)  (70) 

Where  u(x)  is  the  control  vector  and  B(x)  is  the  applications 
matrix.  The  application  matrix,  l>(x),  identifies  the  state 
variables  to  which  u(x)  is  applied.  Like  the  uncontrolled 
constant  coefficient  system,  equation  (70)  cun  be  transformed 
into  a  set  of  uncoupled  differential  equations. 

Defining  modal  variables.  h(x),  su^h  that 

x(x)  =  I-  ( x  )n(x) 

Equation  (70)  can  be  rewritten 

1 '  (x)nix)  +  I- (x ) n  '  (x )  =  A (x )  f  (x ) n  (t )  + 

B(x  )u(x) 

5  4 


.  *  A 


(71) 


S imp 1 i Tying  this  result,  a  set  of  differential  equations  for 
the  modal  variables  is  found, 


n  '  U )  -  f  1  (t)  [A(t)  i:(t)  -  i:,(t)]n(T)  + 
F  1  (  r  )B  ft  )u  (t  J 


(71) 


To  produce  n  uncoupled  differential  equations  before  control 
is  applied,  the  term 

f'1  (t)[A(t)F(t)  -  F'(t)] 


must  be  diagonal.  Recalling  liquation  (59),  Floquet  theory 
provides  the  result, 

F'  (t)  =  A  ( x )  F  ( t )  -  F  (t  )  J 
or  by  rearranging, 

J  =  F'\t)[A(t)1-'(t)  -  1; '  (x )  j  (75) 

Tli  us ,  by  substituting  this  result  of  Floquet  theory  into 
liquation  (72)  the  uncontrolled  decoupled  modal  differential 
equations  are  found, 

h'(t)  =  J’i  (t  )  +  F  1  ft )  B(t  j  u  (t )  (74) 


where  J,  from  Chapter  111,  is  a  Jordan  form  matrix  composed 
of  the  n  characteristic  exponents. 


Mod  e -  font  ro 1 1  a  hi  1 i t  y  Ma  t  r  i  x 

The  mode-controllability  matrix  of  liquation  (74) 
consists  of  two  elements:  the  F  1(t)  and  the  B(t)  matrices. 
Defined  as  the  product, 

S  ;> 


R(t)  ~  I:  (t)R(t) 


(75) 


the  mode-control  lability  matrix  identifies  the  modes  of  the 
system  that  are  controllable  (Ref.  13:  07).  The  matrix  B(t) 
identifies  the  physical  coordinates  to  which  control  is 
applied.  As  such,  it  is,  to  some  extent,  available  as  an 
option  to  the  control  engineer.  However,  the  control  must 
be  applied  to  the  system  in  a  scheme  that  is  physically 
realizable . 

When  the  attitude  of  a  satellite  is  being  controlled, 
the  control  is  applied  to  the  angular  velocity  of  the  body 
about  the  three  axes.  Generally,  this  control  is  applied 
via  gas  jets  or  small  thrusters.  For  t lie  satellite  cases 
discussed  in  Chapter  II,  a  suitable  1>(t)  matrix  to  accom¬ 
plish  this  control,  with  a  scalar  control  law,  is: 


"  0 
!  o 


B(t)=B= 


(  "  b ) 


i  i 


i  J 


Thus,  as  the  control  is  applied,  the  angular  velocity  of  0 

i 

and  6  change,  resulting  in  a  corresponding  change  in 
a n g u  1  a  r  G  i  s  p  1  a c  er.ont . 

The  remaining  component  of  the  mode - contro  1 1  ah i  1  i  t y 
matrix  is  F  1  (x)  ,  which  from  Chapter  III  is  a  periodic 
matrix.  This  matrix  may  be  found  by  integrating  liquation 


(59)  to  obtain  F(x)  and  then  inverting  the  result.  Since 


r 


F  ( T )  must  be  inverted  at  every  instance  control  is  applied, 
this  would  require  extensive  computat  ion  time  or  F  1  i  T )  may 
be  calculated  directly  by  noting: 

i'*1  (T)F(t)  =  L  1J 

Differentiating  this  with  respect  to  x  , 

1;  '(T)';lT)  +  f  =  0 

Substituting  liquation  (59)  into  this  result, 

1;  1  ’  C t )  =  F  1  (  t  )  [  l;  ( x  )  J  -  A  ( x )  F  (t  )  ]  1;  *  1  (t ) 

and  finally,  by  simplifying, 

F'1 ' (t)  =  JF'1 (t)  -  F " 1 (t ) A (t )  (77) 

Integrating  this  matrix  differential  with  the  initial  con¬ 
dition  F  1 (0)  over  one  period  finds  F  5(x)  for  all  x.  Kith 
F"1  (x)  and  B,  the  mode  -  cent ro 1 1 ab  i  1  i  t y  matrix,  is  defined, 
which  since  F  1 (t )  is  periodic,  is  periodic.  Finally,  by 
defining  the  control  Law  u(x),  tiie  periodic  system  can  lie 
controlled  and  further  evaluations  made. 

Font  ro 1  Law 

Since  the  purpose  of  this  controller  is  to  stabilize 
the  attitude  motion,  the  controller  attempts  to  drive  the 
states  to  zero.  To  regulate  this  system,  a  feedback  control 
law  consisting  of  a  gain  matrix,  K,  times  the  modal  vector 


i s  used : 


u  =  Kn  (  t  ) 


(.78) 


where 


K  =  [K  k  k  k  ] 

1  2  3  4 


Since  k  is  a  row  matrix,  the  control  vector  reduces  to  a 
scalar  control  law  of  the  form 


u  =  k  ti  (t)  +  k  n  (x)  +  kn  (t)  +  kn  ( r ) 

1  1  2  2  3  3  3  4 


(79) 


With  this  cont  rol  law,  the  general  form  of  the  closed -- loop 
system  is: 
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n ' (t)  = 


or,  by  simplification 


n  l 


'  i  r  )  = 
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R  k 
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K  k 
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R  k 
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wiiicli  is  in  the  form, 

n  '  U)  =  A*  (t)n(t) 
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where  R-  indicates  the  i1'1  component  of  the  moJc-cont  roJ - 
lability  matrix.  Recalling  from  the  discussion  ol'  the  mode- 
controllability  matrix  that  R(-r)  is  periodic,  the  closed- 
loop  system  given  by  liquations  (SI)  is  also  periodic  in  the 
form  of  liquation  (-15).  As  such,  the  asymptotic  stability 
of  this  closed- loop  system  can  be  established  using  lloquet 
theory.  Thus,  the  performance  of  this  controller  can  be 
readily  evaluated  by  finding  the  eigenvalues  of  the  closed- 
loop  monouromy  matrix. 


Control  Law  implementation 

Using  the  control  law  given  by  liquation  (78),  any 
mode  or  number  of  modes  can  be  controlled.  Control  of  a 
mode  is  accomplished  by  having  a  non-:ero  gain  in  the  K 
matrix  at  a  position  corresponding  to  the  mode  requiring 
control.  Thus,  if  the  second  and  fourth  modes  requires 
control,  that  is,  n  (r)  and  q  (r)  have  undesirable  charac- 

2  4 

teristics,  the  gain  matrix  is 


(82) 


This  would  province  the  closed-loop  system  as  follows: 


n  •  (  i  )  = 


R  k 
i  .’ 

X  +  R  k 


R  k 

3  ^ 

R  k 

4  ? 


R  k 

1  4 

R  k 


R  k 

?  4 

X  +R  k 

4  4  I 


( T )  (83) 


By  inspection,  the  characteristic  exponents  of  this  system 


arc  X  and  A  ,  which  remained  unchanged,  and  the  controlled 

1  3 

characteristic  exponents,  A  *  and  X  *.  bxamining  liquations 

2  4 

(,80)  and  (8  1),  as  the  control  is  applied  to  the  system,  the 
n  uncoupled  differential  equations  are  coupled  by  the  addi¬ 
tional  elements  in  the  columns  of  the  Jordan  form  matrix 
corresponding  to  the  controlled  modes.  Two  general  cases 
appear  when  control  is  applied.  The  first  case  occurs  when 
there  is  one  controlled  mode,  and  the  other  is  when  there 
are  multiple  controlled  modes. 

Single  Controlled  Mode.  When  only  one  mode  is  con¬ 
trolled,  the  closed-loop  system  can  always  be  rearranged  so 
that  the  coefficient  matrix  is  in  either  upper  or  lower 
triangular  form.  In  triangular  form,  the  differential  equat¬ 
ion  for  the  controlled  mode  can  be  solved.  With  this  solution, 
tiie  stability  can  be  determined  directly.  in  fact,  the 
location  of  the  controlled  characteristic  exponent  is  a  linear 
function  of  its  uncontrolled  location,  the  appropriate  clement 
of  the  mode-controllability  matrix,  and  the  gain.  All  of  these 
features  are  demonstrated  in  the  following  example. 

Suppose  the  original  system  had  a  unstable  mode 
correspond  ing  to  r.  (t).  Using  liquations  (78)  and  (7p),  the 
control  law  is 

u  =  h  n  ( x  )  (  8  -1 ) 

3  3 

and  the  closed-loop  system  (from  liquation  (81a))  is, 


•ID 


Rearranging  this  set  of  differential  equations, 
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where  the  coefficient  matrix  is 
As  such,  the  solution  for  q  (t  ) 
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factor.  From  liquation  (.85),  the 
n  it)  is 
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in  upper  triangular  form, 
is  found  using  an  integrating 
differential  equation  for 


n ' (t )  -  [  A  +  R  IT  I  k  ;  q  L :  I  =  0  (.87) 

3  3  3  3  ? 

which  is  a  linear  homogeneous  first  order  differential 
equation.  Using  an  integrating  factor  of  the  form 


n,-(x)  =  exp{-/c  [A  +  R  ('. )k  ]  d  7} 

1  o  '  3  3  3 

an  exact  differential  for  q  (t)  can  lie  found  (Ref.  10:  14). 

3 

Multiplying  liquation  (.87)  by  this  integrating  factor 


T 

n  ’  ( t  )  exp  •(  -  f  (  A  +  R  (  7  Ik  i  d  f  }  -  ex  p 

3 

{-  /  (A  +R  (7 )k  ]  dt, }[A  +  R  (x)k  ]q  (t)  =  0 

0  3  3  3  3  3  33 

which,  hy  inspection,  is  the  exact  differential 


(SS) 


|  T 

17  [n  (T)  expi -/  [A  +  R  (C)k  ]  d  £}]  =  0  (Si)) 

U  l  3  O  3  3  3 

Since  liquation  (S9)  is  an  exact  d  i  f  ferential , 

r 

n  (r)  e\p{-/  [A  +  K  (")k  ]  d  r}  =  C  (90) 

3  o  3  3  3 

where  C  is  the  constant  of  integration.  This  constant  can 
be  found  by  setting  x  equal  to  zero: 


r;  (0)  exp  (0)  =  C 

3 

r;  (0)  =  n 

3  3 

0 

Rearranging  liquation  (90),  the  solution  lor  the  controlled 
mode  is  found, 


x 

n3(x)  =  n3  exp  {;  :>3  *  l^mkj  d  r]  (91) 

o 

Recalling  that  the  mode-cont rol lab i 1 i ty  matrix  is  periodic, 
liquation  (91)  can  be  simplified  further  by  representing 
R  ^  ( 1, )  as  an  infinite  Fourier  series  (Ref.  10:  5S7)  : 
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a  + 
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a  cos 
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n  =  1 


bn  sin 


(92) 


where  the  coefficients  are  given  by  the  Filler  formulas: 
ac  =  1  .•  R3(Od  t 

aa  =  '  I  R 3  ( b )  cos  dr.irxl-T/2,  T/2)J 

and 

b  =  /  R  ( f )  sin  dC 

n  i  3  i 


4  2 


Substituting  liquation  ('JJj  into  liquation  { 9  1  J  one  obtains 


(-)  =  n  exp-.  (\  +  a  k  )t  }  exp xf  R 

3  5  3  0  3  O  3 

0  P  (95 j 

d  m 

3 

where  R  represents  the  periodic  portion  of  R  .  Since  the 

3  p  3 

asymptotic  stability  of  the  closed-loop  system  depends  on 
the  eigenvalues,  or  characteristic  multiples,  of  the  mono- 
dromy  matrix  of  liquation  (S7),  tiie  eigenvalue  associated  with 
the  controlled  node  is  given  by  liquation  (95)  evaluated  at 
t  =  r ,  with  r-;  =1.  evaluating  liquation  (95), 


r,  (i)  =  exp{  f.\  +  a  k  )?}exp{;'  R 

3  3  0  3  0  3 


( •  )  k  d  X  • 


where 


o  3  . 


C,  )k  dr  =  (i 


Simplifying  this  result, 


t i  (?)  =  exp{ (A  +  a  k  )?} 

3  3  0  3 


or,  the  controlled  system  characteristic  exponent  is  given 
by  liquation  (54) 


A  *  =  TT.-expi  (A  +  a  k  )?}  =  A  +  a  k  (< 

3  1  3  C  3  3  0  3 

Thus,  the  closed-loop  characteristic  exponent  is  found  by 
this  linear  equation,  and  with  the  gain  value  such  that 


A  +  a  k  <0 

3  0  3 

stability  is  assured. 


From  F.quation  (97),  control  oi‘  a  single  mode  is  predi¬ 


cated  on  the  constant,  i.e.  zero  frequency ,  portion  of  R  (t) 
being  non-zero.  IF  the  constant  portion  is  equal  to  zero, 
Stability  can  still  be  achieved  by  allowing  the  gain  to  be 


periodic  of  the  form, 

,  .  .  ,  „  .  2iir  i 

k  ( X  )  -  k  C  0  S  —  — 

3  3  1 

(98a) 

or 

k  ( T  )  =  k  sill 

(98b) 

3  3  1 


where  i  is  the  index  of  some  non-zero  Fourier  coefficient  of 
R f  v ) .  Using  this  non-constant  gain  along  with  the  trigonometric 
i dent  it  ites  , 

sin:T  =  hi  1 -cos  2x)  (99a) 

c  o  s  -  t  =  ( 1  +  e  o  s  Z  t  )  (99b) 

a  non- zero  constant  portion  of  R  (x)  is  generated.  This 
procedure  is  demonstrated  as  follows,  using  liquation  (98a), 

k  (r)  =  k  cos  —  —  j 

3  3 

Multiplying  this  result  by  the  appropriate  clement  of  the 
mode-controllability  matrix,  one  finds, 

R  (t)  k  (x)  =  a  k  cos  j  +  (a;cos  r- i ) 

3  3  03  1  1-11  i 

i  i  j 

(k  cos  )  +  ■  :=  !  ( b  |  s  i n  A.- i )  (  k  cos  1  (100) 

*  ^1. 


+  a . k  cos 


With  the  last  term  ot  this  result  and  liquation  (99b),  a 
constant  portion  is  generated: 


a  .  k 


I  2  -  ”  T  •  I  3/i  4  r  T  • 

.  k  cos'  — -J  =  -^—(1  +  cos  -v-  j) 


J  3 


or  the  constant  portion  is 
a  . 

,  -  J 


With  the  result,  the  controlled  mode  is  always  controllable 
provided  R  (tl  is  not  zero  for  all  t.  Otherwise,  if  R  (x)  = 

3  ? 

0,  the  third  mode  is  uncontrollable.  This  is  the  result  that 
the  inode-controllability  matrix  indicates. 

Thus,  when  there  is  one  controlled  mode,  the  solu¬ 
tion  to  the  closed-loop  periodic  differential  equation  can 
be  solved  using  an  integrating  factor.  Using  this  result, 
the  stability  of  the  closed-loop  controlled  mode  can  be 
determined  directly  by  using  liquation  (96).  Since  modal 
control  only  controls  the  desired  mode(s),  the  stability 
of  the  other  modes  is  unaffected  by  the  control,  and  the 
stability  of  the  entire  system  is  assured. 

Multiple  Controlled  Modes.  When  there  are  multiple 
controlled  modes,  the  stability  is  not  as  readily  attainable  as 
it  was  in  the  previous  case.  The  reason  for  this  is  that  the 
closed- loop  coefficient  matrix  can  not  be  rearranged  to  form 
a  triangular  form,  either  upper  or  lower,  matrix.  Since  the 
matrix  is  not  in  triangular  form,  the  solution  to  the  differ¬ 
ential  equations  for  the  controlled  modes  arc  not  readily 


obtained . 


As  an  example,  suppose  the  original  system  had  two 


unstable  modes,  n  (t)  and  n  (x).  From  liquation  (Sla),  the 

2  3 

closed -loop  system  is  represented  by 


A  '  (.  t  J 


R  k 
1  2 

X  +R  k 


R  k 

3  2 


R  k 

1  3 

R  k 


X  +  R  k 

3  3  3 


R  k 


R  k 


(x)  (101) 


Unlike  the  case  when  only  one  mode  is  controlled,  the 
coefficient  matrix  of  this  equation  can  not  be  rearranged 
into  a  triangular  form.  Therefore,  to  find  the  closed-loop 
characteristic  exponents  a  Tloquet  theory  analysis  using 
the  concepts  presented  in  Chapter  ill  must  be  completed. 

By  inspection  of  liquation  (101),  the  closed-loop 
characteristic  exponent  arc  X  ,  X  and  the  two  found  from 

i  *t 

the  monodromy  matrix  calculated  by  integrating  the  sub¬ 
matrix  , 


4'(x) 


R  k  +X 
2  2  2 


<Hx) 


R  k 


3  2 


X  +R  k 

3  3  3 


There  are . essential ly  two  methods  of  setting  the  gains  to 


stabilize  this  s  vs  tort.  Tin 


e  t  no  a  c  on  s  i  s  t  s  o  t  a 


systematic  search  using  various  gain  settings  until  the 
desired  characteristic  exponents  are  found.  The  second 
method  entails  successive  diagonal icut ion.  In  this  method, 


one  gain  is  set  resulting  in  a  controlled  system.  Then 


this  controlled  system  is  rediagonal i zed  to  form  a  new 
Jordan  form  matrix  and  another  gain  is  then  set.  This  process 
is  repeated  until  all  modes  are  stabilized. 

This  completes  the  development  of  the  control  law 
and  the  mode-controllability  matrix.  This  control  law 
presented  in  this  chapter  is  capable  of  controlling  any 
number  of  modes.  in  the  following  chapter,  example  test 
cases  from  both  the  symmetric  and  unsymmetric  satellites 
are  used  to  demonstrate  the  actual  usage  of  this  control 
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RESULTS 

In  Chapter  II  of  this  study,  the  dcvcl  opment  of 
linearized  equations  for  the  attitude  motion  of  two  satel¬ 
lites  was  presented.  In  both  cases,  the  stability  of  these 
periodic  systems  was  found  to  be  dependent  on  various 
parameters,  including  the  orbital  eccentricity,  moment  of 
inertia  ratios,  and  spin  rates.  Obviously,  by  varying 
these  parameters,  an  infinite  number  of  test  cases  are 
found.  Therefore,  five  example  test  cases  were  selected 
from  these  two  satellite  cases.  Using  these  five  test  cases, 
the  theoretical  developments  of  this  study  are  verified  with 
numerical  results.  Upon  selecting  the  live  test  cases,  a 
stability  analysis,  using  l'loquot  theory  and  the  results 
of  Chapter  III,  is  conducted.  Kith  this  stability  analysis, 
two  of  the  test  cases  arc  digitally  simulated  to  verify 
the  results.  Next,  three  of  the  test  cases  are  controlled 
using  the  developments  of  Chapter  IV  and  the  closed-loop 
linear  system  is  simulated.  finally,  numerical  problems 
which  appeared  several  times  are  presented  using  one  of  the 
test  cases. 

Test  Cases 

from  the  infinite  number  of  combinations  of  satel¬ 
lite  parameters,  five  test  cases  arc  chosen  for  use  during 

•  IS 


the  presentation  of  the  numerical  results.  Of  these  live, 
two  arc  symmetrical  satellites  in  elliptical  orbits,  and 
the  remaining  three  are  unsymmetr i cal  satellites  in  circular 
orbits.  Table  111  is  a  summary  of  these  five  test  cases  and 
their  associated  parametric  values.  The  test  cases  numbers  arc 
used  throughout  this  chapter  to  identify  the  individual  cases. 

Table  111 

Parametric  Values  of  the  Five  Test  Cases 


Test  Case 

e 

a 

k 

i 

k 

2 

1 

0.5 

-1.0 

1  .  0 

— 

? 

0.5 

-2.1) 

.  (> 

— 

3 

0.0 

-1.0 

"> 

- .  0 

4 

0.0 

I  .  0 

- .  0 

.  0 

5 

0.0 

1  .  0 

-1.0 

-  .  3 

Stability  Analysis 

The  first  step  in  the  control  of  an  arbitrary  satel¬ 
lite  design  is  to  conduct  a  stability  anaylsis  of  the  open 
loop  system.  When  the  linear  system  has  constant  coeffi¬ 
cients,  the  stability  is  easily  determined  by  the  eigenvalues 
of  the  coefficient  matrix.  However,  when  the  system  is 
periodic,  Floquet  theory  is  used  to  ascertain  the  asymptotic 
stability.  Using  the  theory  of  Chapter  III,  the  stability 


of  the  five  satellite  eases  is  determined.  Tabic  IV  contains 
the  four  characteristic  exponents  found  using  floquet  theory 
for  each  test  case.  from  Chapter  111,  asymptotic  stability 
is  assured  if  the  real  part  of  all  characteristic  exponents 
is  less  than  zero. 


Table  IV 

Characteristic  exponents  of  the  Five  Test  Cases 


Case 

X 

l 

X 

3 

X 

u 

1 

.0718 

-.0718 

0 . 0 

0.0 

+  .  3728  i 

- . 37  2  S  i 

T 

u 

0.0 

0.0 

0.0 

0.0 

+  . 1  2  00  i 

- . 1 20Oi 

' . 4  00  5 i 

- .4005 i 

3 

.  12  58 

- . 1258 

0 . 0 

0.0 

+.18101 

-  .  1 8 1 9  i 

4 

.  3  5  5  T 

.  3  5  5  7 

-.5557 

- . 3357 

+  .  1  7  7  o  i 

-  .  1  7  7  0  i 

+  .  1  7  7  0  i 

-  .  1  7  7  0  i 

3 

.  8  2  o  5 

-  .  8  2  o  5 

.2128 

.2128 

When  the  character i st i c  exponents  have  a  real  part  of  aero, 
the  system  exhibits  neutral  stability.  Of  the  five  test 
cases,  all  are  unstable  except  Case  2.  Case  2  is  stable,  and 
it  exhibits  the  characteristics  of  stable  gravity  gradient 
satellites  mentioned  in  Chapter  1.  With  the  characteristic 
exponents  for  Case  2,  a  stable  satellite  has  eero  damping  and 
low  natural  frequencies  of  oscillation.  Thus,  the  four  modes 


of  this  stable  case  arc  pure  oscillatory:  very  undesirable 


attitude  motion  for  a  precise  satellite  system.  Cases 

1  and  5  have  one  unstable  mode,  and  Cases  4  and  5  have  two 
unstable  modes.  Cases  1  and  5  both  have  one  unstable  real 
characteristic  exponent,  Case  4  has  two  unstable  complex 
conjugate  characteristic  exponents,  and  Case  5  lias  two  un¬ 
stable  real  characteristic  exponents. 

To  validate  the  results  of  the  stability  anaylsis, 
the  uncontrolled  linear  systems  of  Cases  2  and  5  were  simu¬ 
lated  over  fifteen  periods.  At  first,  the  uncontrolled 
system  of  Case  2  was  simulated  using  physical  coordinates 
as  the  basis  variables.  The  result  are  presented  in  figures 

2  through  5.  As  expected,  no  specific  trends  on  asymptotic 
stabiltiy  are  ascertained  from  these  figures.  The  reason 
for  the  simulation  results  appearing  as  they  do  is  that  the 
state  variables  are  actually  linear  combinations  of  the  four 
modes.  The  characteristic  exponents  represent  the  response 
of  the  individual  modes,  not  the  physical  state  variables, 
from  liquation  (71), 

Mr)  =  1:  {  t  )  n  1  t  ) 

and  the  state  variable  responses  are  linear  combinations  of 
the  responses  of  the  individual  modes.  Therefore,  to  gain 
more  insight  into  the  responses  of  the  various  test  cases 
when  simulated,  all  simulations  results  are  presented  in 
modal  variables. 

The  initial  conditions  used  throughout  this  study. 


a  I 


for  simulation  purposes  are: 


Figure  2.  Uncon t rol 2 cd  x  ft)  Response  for  Cose 


figure  5.  Uncontrolled  x  (x)  Response  for  Case 


-urn  = 

This  initial  conditions  vector  was  chosen  to  assure  that 
all  modes  are  excited.  Also,  all  simulation  results  are  over 
a  total  of  fifteen  periods.  This  period  length  was  assumed 
to  be  of  sufficient  length  to  illustrate  the  general  trends 
in  stability  of  the  linearized  equations . 

Simulating  Case  2  in  modal  coordinates,  a  response 
typical  of  a  stable  gravity  gradient  satellite  is  obtained. 
The  responses  of  the  individual  modes  to  the  input  initial 
conditions  are  presented  in  Figures  6  through  9.  As  ex¬ 
pected,  the  first  and  second  modes  represent  the  extremely 
low  frequency  oscillations  and  the  third  and  fourth  modes 
represent  the  slightly  higher  frequency  oscillations.  Also, 
from  Figures  8  and  9,  no  damping  of  the  oscillatory  motion 
is  observed. 

To  verify  the  instability  caused  by  positive  real 
parts  of  the  characteristic  exponents,  Case  5  was  simulated. 
From  Table  IV,  the  first  mode  is  unstable.  Figures  10 
through  3  depict  the  responses  of  the  four  modes  of  Case 
3.  As  expected,  the  first  mode,  the  mode  with  positive  real 
characteristic  exponent,  is  definitely  unstable.  The  stable 
real  mode,  the  second  mode,  demonstrates  the  typical  response 
for  a  mode  with  a  negative  real  portion.  Finally,  the  third 

.">() 


Figure  6.  Uncontrolled  I:irst  Mode  Response  for  Gas 


:igure  8.  Uncontrolled  Third  Mode  Response  for  Case 


controlled  Fourth  Mode  Response  for  Case 


Uncontrolled  First  Mode  Response  for  Case 


Uncontrolled  Second  Mode  Response  for  Cose 


Uncontrolled  Fourth  Mode  Response  for  Case 


and  fourth  modes  represent  the  typical  low  frequency,  low 
damping  oscillatory  response  of  neutraLly  stable  complex 
conjugate  modes. 


Thus,  using  Floquct  theory  as  presented  in  Chapter 
III,  the  asymptotic  stability  of  the  five  test  cases  is 
determined.  By  simulating  the  uncontrolled  system,  the 
stability  of  the  uncontrolled  systems  is  verified.  In  the 
following  section,  three  of  the  unstable  satellite  cases  arc 
stabilized  using  the  modal  control  theory  of  Chapter  IV. 

Satellite  Stabilization  and  Verification 


In  this  section,  three  of  the  four  unstable  satellite 
test  cases  are  stabilized  using  modal  control  techniques  and 
the  control  law  developed  in  Chapter  IV.  These  three  cases, 
cases  1,  4,  and  5,  represent  three  different  instability 
configurations.  Case  1  lias  one  unstable  real  characteristic 
exponent.  Cases  4  and  5  both  have  two  unstable  modes.  Case 
4  has  two  unstable  complex  conjugate  characteristic  exponents 
and  Case  5  has  two  unstable  real  characteristic  exponents. 

Single  Mode  Stabilization.  To  control  Case  1,  a 
controller  must  be  utilized  to  control  only  the  first  mode, 
therefore,  tiie  results  from  the  Single  Con  t  ro  1  1  ed  Mode  sec¬ 
tion  of  the  previous  chapter  are  used.  Since  there  is  one 
controlled  rode,  the  control  law  is 

'll  r  )  =  hr)  It) 


From  Chapter  IV,  liquation  (96)  states  that  the  closed  -  loop 
characteristic  exponent  is 

X*  =  X  +  a  k 
1  1  0  1 

For  this  case,  the  uncontrolled  system  characteristic 
exponent  is  0.6718  and  the  zero  frequency  portion  of  the  first 
element  of  the  mode-controllability  matrix  is  0.6198.  There¬ 
fore,  the  gain  required  to  move  the  characteristic  exponent 
to  zero  is 

0.0  =  .6178  +  .6198k 

i 

o  r 

k  =  -1.084 

i 

For  the  actual  control  of  this  system,  the  gain  was  set  at 
-2.0,  which  results  in  a  controlled  characteristic  exponent 
of  -.5677.  Table  V  contains  the  controlled  system  charact¬ 
eristic  exponents  for  Case  1  based  upon  the  numerical  Floquct 
analysis  described  in  Chapter  III. 

Table  V 

Case  1  Characteristic  exponents  with 

K=  (-2.0,  0.0,  0.0,  0.0) 

Uncontrolled  Cc  lied 


.  o  7  1  8 
-.0718 

0.0  +  .  3 7  2  8 i 
0.0  -  .57281 


-  .  o  7  1  8 

0.0  +  .  5  7  2  8  i 
0.0  -  .  5 7  2  8  i 


as  expected,  the  simulation  results,  Figures  14  through  17, 
show  a  satellite  with  two  exponentially  decaying  modes  and  two 
oscillating  modes.  lixamining  these  responses,  one  notices  a 
phenomenon  that  did  not  appear  in  the  uncontrolled  simulations, 
liach  mode  contains  a  transient  response  that  decays  rapidly. 
This  transient  response  is  due  to  the  controller  modifying 
the  eigenvector  associated  with  the  controlled  mode  of  the 
closed-loop  monodromy  matrix.  This  new  eigenvector  causes 
the  initial  transient  to  be  displayed  in  all  modal  responses. 
With  the  addition  of  the  controller,  the  characteristic 
exponents  associated  with  the  second,  third,  and  fourth 
modes  remained  unchanged.  Only  the  characteristic  exponent 
for  the  first  mode  was  altered.  These  results  verify  that 
modal  control  techniques  only  .affect  the  character ist  ic 
exponents  of  the  controlled  modes.  Also,  the  numerical 
results  of  this  example  validate  the  theory  concerning  the 
control  of  one  mode  in  Chapter  IV. 

Multiple  Controlled  Mode  Stabilisation .  Roth  cases 
4  and  5  have  multiple  unstable  modes.  From  Chapter  IV,  we 
found  that  when  there  are  multiple  controlled  modes,  the 
gains  required  for  stabilization  cannot  be  found  directly 
by  using  a  linear  formula  like  the  previous  example.  This 
problem  with  gain  selection  was  due  to  t  he  cent ro 1  led - syst cm 
coefficient  matrix  not  being  in  a  triangular  form.  There¬ 
fore,  the  solution  to  this  type  of  problem  is  accomplished 
using  a  systematic  search  using  various  gain  settings. 

This  technique  was  used  to  stabilize  both  Case  4  and  Case  5. 
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Svstem  Second  Mode  P.cs 


Controlled  System  Third  Mode  Response  for  Case 


rolled  Systcr.  fourth  Mode  Response  for  Cnsc 


From  Table  IV,  the  unstable  modes  of  Case  4  are  the 
lirst  and  second  modes.  Therefore,  the  general  control  law 
necessary  to  control  this  case  is 

’‘It )  =  h  n  (t )  +  k  r,  (t) 

II  2  2 

Through  the  use  of  a  systematic  search  procedure,  the  values 

of  the  two  gains  which  provided  adequate  control  and  stable 

system  characteristic  exponents  are  k  =  4.0  and  k  =  -10.0. 

1  2 

With  these  two  gains,  the  two  unstable  characteristic  ex¬ 
ponents  were  moved  to  the  left  side  of  the  complex  plane. 

The  results  oi  the  controlled-system  monodromy  matrix  cal¬ 
culations  are  presented  in  Table  VI.  From  these  results, 
once  again  modal  control  techniques  controlled  only  the 
modes  requiring  control. 


Table  VI 

Case  4  Characteristic  exponents  with 


K= 

(4 

.0,  -10.0, 

0.0,  0.0) 

Uncontrolled 

Control  led 

x_ 

X 

.  555. 

+ 

.  1  7  7  6  i 

-1.011  + 

.  0 1 4  0  i 

X 

2 

.  5557 

+ 

.  17  761 

-1.011  - 

.  0 1 4  b  i 

A 

3 

-.5557 

+ 

.  1 7  7  6  i 

-  .5557  + 

.  1  7  7  0  i 

A 

-.5557 

- 

.  1 7  7  6  i 

-.5557- 

.  1  7  7  6  i 

These  gain  settings  were  used  in  a  digital  simulation 
of  the  modal  variables.  The  results  of  the  simulation  are 


illustrated  in  Figures  18  through  21.  In  this  case,  the 
oscillations  in  all  modes  damp  out  within  three  periods. 

This  is  due  to  the  location,  and  therefore  the  existence  of 
other  than  negligible  damping,  in  the  complex  plane,  of 
the  controlled  system  characteristic  exponents.  Once  again, 
the  initial  transients  in  all  modes  due  to  the  new  eigen¬ 
vectors  of  the  controlled  system  are  apparent  in  these 
figures.  Since  all  modes  of  this  system  damp  out  quickly 
and  there  are  no  highly  oscillatory  modal  responses,  this 
system  presents  an  excellent  system  for  simulation  in 
physical  variables.  Therefore,  with  the  initial  conditions, 


x(0)  = 


.  05326  | 
.  0  5326 
0.0  I 

I 

0.0  I 


and  the  control  law  given  by 


u(t)  =  KF  1  (t)x(t) 


this  controlled  system  was  simulated  for  fifteen  periods 
using  the  physical  variables.  This  simulation  verifys  that 
the  linearized  physical  system  is  controlled  with  no  un¬ 
bounded  growtli  in  time  of  any  of  the  state  variables.  Figures 
22  through  25  illustrate  the  results  of  this  simulation. 

The  controller  performed  as  expected,  all  oscillations  of 
the  state  variables  stopped  after  approximately  three 
periods . 


Figure  IS.  Controlled  System  First  Mode  Response  for  Case 


Figure  21.  Controlled  Systen  Fourth  Mode  Response  For  Case 


Controlled  x  ft)  Response  for  Case 


Figure  25.  Controlled  x  (t)  Response  for  Case 
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Controlled  x  (x)  Response 


SI 


The  final  test  ease  to  he  controlled,  Case  5,  lias 
two  unstable  real  characteristic  exponents.  Thus,  since 
there  are  two  modes  requiring  control,  the  same  techniques 
used  in  the  control  of  Case  4  are  used.  From  Table  IV,  the 
unstable  inodes  are  the  first  and  third  modes.  Therefore, 
the  scalar  control  law  for  the  control  of  this  system  is 

u(t)  =  k  n  (T)  +  h  n  (t) 

11  3  3 

Using  a  systematic  search  technique,  the  gains  that  provided 
adequate  control  for  this  satellite  are  k  =  -1.4  and  k 

1  3 

-.01.  Table  VII  contains  the  results  of  the  closed- loop  ana¬ 
lysis  using  Floquet  theory. 


Table  VII 

Case  5  Characteristic  exponents  with 
k=  (-1.4,  0.0,  -.01  ,  0.0) 


Uncontrolled 

Cont rol 1 cd 

X 

1 

.  8265 

- .  058  2  +  .  097  1  i 

A 

2 

- .8265 

-.8265 

A 

3 

.2128 

-.  0582  -  .  0  9  7 1 i 

A 

4 

-  .  2128 

- . 21 28 

From  these  results,  the  control  law  provides  adequate 
control  for  the  unstable  modes  to  bring  them  under  control. 
Figures  26  through  29  present  the  results  of  simulating  the 


Controlled  System  First  Mode  Response  for  Car 


27.  Controlled  System  Second  Mode  Response  for  Case 
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Controlled  System  Third  Mode  Response  for  Cose 


Controlled  System  Fourth  Mode  Response  for  Case 


controlled  system  over  fifteen  periods.  In  this  case,  all  modes 
were  controlled  and  eventually  reduced  to  zero.  The  previously 


mentioned  transient  responses  are  more  pronounced  in  this  case 
than  in  the  other  two  cases.  These  transient  response  char¬ 
acteristics  are  due  to  the  eigenvectors  of  the  closed-loop 
system  monodromy  matrix.  During  this  study,  no  attempt  was 
made  to  specify  the  desired  transient  response  characteristics 
of  the  controlled  systems. 

With  these  numerical  and  graphical  results,  the 
implementation  and  verification  of  the  control  law  and 
associated  theoretical  developments  arc  complete.  As 
presented  in  the  examples,  modal  control  provided  the  required 
control  over  only  the  controlled  modes  and  only  the  initial 
transients  of  the  other  modes  were  affected.  Also,  the  general¬ 
izations  made  concerning  stable  gravity  gradient  satellites 
in  Chapter  I  were  demonstrated  in  the  various  examples. 

During  the  completion  of  these  investigations,  numerical 
problems  associated  with  the  computations  of  the  F  1  (t ) 
matrix  were  encountered.  In  the  fo; lowing  section,  these 
numerical  problems  and  their  effects  on  the  controlled  system 
are  addressed. 

Numerical  Problems  Associated  with  F  1 (i 1  Computations 

One  of  the  basic  elements  in  the  closed- loop  system, 
from  Equation  (,74),  is  the  F  1  (t )  matrix.  This  matrix,  which 
appeared  during  the  d iagonal i zat ion  and  subsequent  control 
processes,  is  periodic  with  period  T.  As  such,  by  finding 
F  1 (t )  over  one  period  is  equivalent  to  finding  this  matrix 

8  7 

I 
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over  all  x.  Therefore,  to  case  computational  burden,  the 
method  of  computing  the  i:  l(x)  matrix  for  the  closed-loop 
stability  analysis  and  subsequent  digital  simulations  was 
to  compute  F  1(r)  for  m  evenly  spaced  points  during  one 
period.  Using  these  in  F  '  (x)  matrices,  a  harmonic  analysis 
was  completed  on  all  sixteen  elements  (Ref.  2:  108-109). 

The  results  of  this  harmonic  analysis  were  ^  +  1  cosine 
Fourier  coefficients  and  ^  -  1  sine  Fourier  coefficients. 

With  these  coefficients,  the  F  ’(x)  matrix  can  be  computed 
at  any  x.  Using  this  technique,  the  closed- loop  system,  in 
either  modal  or  state  variables,  is  integrated.  Using 
Case  3,  it  is  shown  how  this  computational  technique  does  not 
insure  the  precise  computations  of  the  F_1(x)  matrix  required 
to  satisfy  equation  (77).  From  Table  IV,  Case  3  has  one  un¬ 
stable  characteristic  exponent.  Thus,  since  there  is  one  un¬ 
stable  mode,  the  controlled  characteristic  exponent,  A*,  is 
found  by  using  Equation  (96).  For  this  case,  the  zero  fre¬ 
quency  portion  of  the  first  clement  of  the  mode-controll¬ 
ability  matrix,  a  ,  is  -0.2187.  Using  Equation  (96),  the 

0 

theoretical  value  for  A*,  is  found,  and  by  finding  the 
characteristic  exponents  of  the  closed- loop  system  monodromy 
matrix,  the  numerically  derived  value  of  A*  is  found.  Com¬ 
paring  these  two  values,  the  effects  of  these  numerical 
problems  associated  with  the  F  J(x)  can  be  ascertained. 

These  computations  were  carried  out  for  five  gain  settings, 
the  results  arc  tabulated  in  Table  VIII.  From  Table  VIII, 
one  notices  that  as  the  gain  is  increased,  the  theoretically 


Theoretical  and  Computed  Values  of  the  Characteristic 


and  numerically  computed  characteristic  exponents  diverge. 

In  fact,  there  should  he  no  variation  in  the  characteristic 
exponents  of  the  modes  which  are  uncontrolled.  However, 
with  a  gain  of  1.0,  there  are  apparent  variations  in  the 
uncontrolled  characteristic  exponents. 

These  variations  are  due  to  the  loss  in  accuracy  of 
1'  1  ( x )  .  To  examine  this  loss  of  accuracy,  the  T  1  (.  t  )  matrix 
at  x  =  it  was  computed  using  the  Fourier  series  for  each 
element  and  by  integrating  liquation  (77)  as  presented  in 
Chapter  IV.  The  results  are  presented  in  Figure  50.  exam¬ 
ining  these  two  matrices,  one  notes  that  there  arc  significant 
differences.  The  accuracy  of  the  Fourier  series  can  be 
improved  by  increasing  the  number  of  function  evaluations 
used  in  the  harmonic  analysis.  An  increased  number  of 
function  evaluations  was  tried,  but  the  resulting  F  1  (x) 
matrix  was  still  too  inaccurate.  Therefore,  these  inaccuracies 
are  due  to  using  the  Fourier  coefficients  to  compute  t  It  c  F  1  (.  x ) 
matrix.  there  are  many  possible  causes  for  these  numerical 
problems,  among  these  are  truncation  errors  and  finite 
wordlength  of  the  computer.  However,  no  definite  cause  for 
these  problems  was  determined. 

This  completes  tire  discussion  of  the  results  com¬ 
piled  during  this  study.  Consisely,  the  results  indicate 
the  ability  to  control  unstable  satellites  under  the 
influence  of  a  gravity  gradient.  However,  numerical 
problems  in  the  computation  of  the  mode-control labil t iy 
matrix  were  encountered. 
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Figure  50.  F  *(x)  Calculated  at  x  =  it 

using  the  Fourier  Coefficients 
and  Analytical  Techniques 


CHAPTER  VI 


CONCLUDING  REMARKS  AND  Rl-COMMliK  DAT  IONS 

The  problem  that  this  study  addresses  is  the  control 
ol  satellite  attitude  motion  '  lie  the  satellite  is  under 
the  influence  of  a  gravity  gradient.  Since  the  linearized 
equations  describing  the  attitude  motion  are  periodic,  the 
controller  must  control  the  location,  in  t lie  complex  plane, 
of  the  characteristic  exponents.  As  seen  in  the  previous 
chapter,  the  modal  control  scheme  and  the  control  law 
developed  during  this  study  provide  adequate  control  over 
the  satellite  attitude  motion.  The  control  law  demonstrated 
that,  in  general,  it  only  controls  the  modes  requiring 
control.  Therefore,  the  modal  control  techniques  provide 
a  viable  control  scheme  for  the  attitude  control  of  un¬ 
stable  satellites  described  by  linear  periodic  systems. 

During  this  study,  it  became  apparent  that  there 
are  two  areas  requiring  additional  work.  first,  a  problem 
was  encountered  numerous  times  with  the  computations  of  the 
mode-controllability  matrix.  The  result  of  this  computation¬ 
al  problem  is  inaccurate  or  fictitous  results  when  the 
computing  the  stability  characteristics  of  the  controlled 
system.  The  second  area  requiring  additional  study  is  in  the 
gain  selection  techniques  for  stabilizing  a  satellite  with 
more  than  one  unstable  mode.  In  this  study,  a  systematic 
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search  procedure  was  used,  hut  in  general,  one  would  like 
to  select  the  gains  directly.  With  additonal  work  in  these 
areas,  the  modal  control  technique  and  the  control  law  pre¬ 
sented  in  this  study  will  provide  an  extremely  flexible 
method  for  controlling  periodic  systems  whose  asymptotic 
stability  characteristics  are  undesirable. 

Recommendations  for  future  Studies 

The  recommendations  for  future  studies  consist  of 
two  levels.  The  first  level  recommendations  are  future  study 
efforts  involving  the  same  techniques  used  during  this  study. 
The  second  level  recommendations  are  for  future  studies  for 
developing  a  control  scheme  for  on-line  applications  in  a 
real  world  environment. 

The  first  level  recommendations  are  ways  to  obtain 
more  accurate  calculations  of  the  mode-controllability 
matrix.  There  is  one  basic  means  of  accomplishing  this. 

This  means  is  by  using  double  precision  arithmetic  F  1 (x) 
calculations.  All  computations  in  this  study  were  completed 
using  single  precision  arithmetic.  By  increaseing  the  pre¬ 
cision,  more  accurate  results  may  be  obtained. 

The  second  level  recommendations  arc  possible  methods 
of  achieving  acceptable  results  in  a  real  world  envoronment. 
The  first  method  is  to  use  a  Kalman  filter  to  provide  state 
estimations.  Using  sensor  models,  which  take  noise  corrupted 
measurements  of  output  variables,  the  Kalman  filter  can 
then  be  used  to  generate  state  estimates.  These  state 


estimates  arc  then  used  as  the  feedback  variables  instead  of 
the  actual  states.  This  technique  provides  a  more  realistic 
application  of  the  control  techniques  presented  in  this 
study.  In  the  real  world,  sensors  would  be  used  to  measure 
the  angular  rates  and  angles.  This  data  would  then  he  used 
to  activate  the  controller.  The  second  method  is  to  design 
an  optimal  controller  using  IaHI  design  techniques.  TOC 
design  techniques  provide  an  optimal  controller  for  a  pro¬ 
blem  described  by  linear  models  with  gradratic  cost  criterion 
and  Gaussian  noise  models.  This  is  a  reasonable  choice 
since  the  linear  periodic  system  models  can  be  assumed  to 
1)0  corrupted  by  noise.  This  noise  can  be  due  to  such  items 
as  atmospheric  drag,  gravity  anomalies,  solar  wind,  etc. 

Under  these  assumptions,  an  optimal  stochastic  regulator 
can  be  designed. 
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